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I nt  rcduct ion 


I  . 


a.  Scope.  The  purpose  of  this  report  is  to  present  an  engineering 
method  of  evaluating  sky-wave  field  intensities  in  tne  high  freai-ency 
band  *or  any  t ransm i ss i on  distance  cp  to  about  Ip, 000  kilometers. 
Field  intensity  recordings  were  statistically  analyzed  in  order  to  study 
daytime  absorption,  day-to-day  variations,  and  absolute  magnitudes. 
Based  upon  the  analysis,  t  proposeo  prediction  method  was  obtained  which 
should  give  relatively  uccuraie  sky-wave  lield  intensities  for  the  fre¬ 
quencies  and  distances  considered.  Succeeding  sections  show  the  method 
of  analysis  of  absorption,  the  method  of  analysis  of  absolute  magni¬ 
tudes,  the  miethod  of  analysis  of  day-to-day  variations,  the  prooosed 
fiald  intensity  prediction  method,  and  summarizing  conclusions.  All 
necessary  equations  used  in  the  prediction  netnod  are  represented  in 
nomograms. 

h.  DATA.  The  data  analyzed  consists  of -monthly  median  receiver  input 
voltages  at  one  hour  time  intervals  throughout  the  day  and  mass  plots 
shewing  hourly  median  receiver  input  voltages  i a  2  decibel  class  inter¬ 
vals  for  the  days  of  the  month.  The  circuits  analyzed  arc  listea  in 
Table  I.  This  table  snows  the  perioo  for  which  the  data  was  available, 
the  distance  between  the  transmitting  and  receiving  station,  the  wave 
frequency,  and  type  of  data.  Tnese  records  were  continuous  over  a  peri¬ 
od  of  at  least  as  long  as  one  year.  The  circuit  distances  rar.ge  from 
bp  i o  15,000  kilometers  and  the  frequencies  from  0.7  to  20  megacycles 
por  second.  Eighty-three  circuit  years  of  monthly  median  hourly  median 
receiver  input  voltages  were  used  in  the  daytime  absorption  analysis. 
Data  fo'  S6  circuits,  including  information  about  transmitter  output 
powers,  transmitting  antenna  ga:ns,  and  receiving  system  input  charac¬ 
teristics  were  used  i n  the  analys i s o f absol ute  magnitudes.  The  receiving 
system  installations  «eie  so  designed  as  to  make  estimations  of  absolute 
magnitudes  possible.  Forty-three  years  of  massolot  data  was  used  in  the 
analysis  of  day-to-day  variations  of  hourly  median  field  intensities. 
Other  data  not  shown  in  the  table  covering  a  period  from  1037  to  1944. 
was  used  in  studying  the  relationship  between  daytime  absorption  and  !2 
month  running  average  sunspot  numter. 

c.  The  Method  or  Analysis.  This  analysis  isbased  upon  tne  assumption 
that  at  hiqh  frequencies  sky-wave  field  inxensities  car.  oe  computed  by 
geometric  optical  methods.  Various  authors  show  the  applicability  of 
.ay  theory  at  these  frequencies.1  Sky  wave  propagation  is  a  geometric 
rather  than  a  diffraction  problem  since  each  place  on  earth  at  distances 
qreater  than  the  skip  distance  is  accessable  by  georetric  t  raj  ect  or  i  ef 
\  method  such  as  a  residue  seiies  would  be  mo'e  apolicable  if  a  large 
number  of  geometric  rays  existed  all  having  the  same  order  of  intensity. 
It  has  been  shown  by  L.  E.  Beqhian  that,  for  high  frequencies,  only  a 
few  rays  are  required  and  that  under  conditions  of  absorption  the  first 
ray  approximates  the  field  strength  obtained  by  assuming  ail  possible 
rays. ^  When  several  ray  paths  exist  having  the"  same  order  of  magnitudes, 
it  is  necessary  to  evaluate  the  field  by  summing  the  contribution  of  all 
the  rays.  Occasionally  for  long  distances,  the  error  in  assuming  that 
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cne  field  intensity  can  be  represented  ov  one  ray  alone  can  be  relatively 
large  since  ray  intensit ies of  equal  magnitudes  are  possi b!e  along  several 
ray  paths. 

2.  The  Analysis  of  Sky-wave  Absorption 

a.  Gt.  n  t  ft  a  i  .  (I)  In  the  ionosphere,  gas  molecules  are  ionized  into 

electrons  and  positive  :ons  principally  by  the  action  of  ultra-violet 
solar  light.  A  passing  radio  wave  sets  electrons  in  vibratory  motion. 
The  elections  randomaly  collide  with  gas  molecules  and  lose  part  of  the 
energy  of  vibration  during  collisions.  The  absorption  of  a  radio  wave 
per  unit  path  length  depends  on  the  collision  frequency  and  the  number 
of  electrons  per  unit  volume,  both  of  wnich  vary  with  height  above  earth. 

i2)  During  aay'iqiii  hours,  the maj or i tv  of  sky-wave  absorption  at  high 
frequencies  occurs  in  the  lower  E  region.  This  region  is  known  as  the  0 
region  although  it  is  not  certain  that  it  is  physically  distinct  from  the 
c.  The  collision  frequency  is  high  while  the  ion  density  is  low  so  that 
a  radio  wave  is  absorbed  but  not  greatly  refracted.  This  type  of  absorp¬ 
tion  is  called  non-devlatl ve  absolution.  The  geometric  ray  theory 
approximation  is  app  I  i  cab  I  e  fo  r  f  requenc  i  es  nottooclose  to  the  collision 
frequency,  the  ionosphere  approximating  a  part  i  a  I  ly  ref  I  ect  i  ng  dielectric. 

< 3 >  During  the  night  hours,  the  ionosphere  acts  as  an  almost  totally 
reflecting  dielectric  at  high  frequencies.  The  effects  of  the  passage 
of  radio-waves  through  the  0  and  E  layers  are  neg I i g i bl e,  since  the  ioni¬ 
zation  of  these  layers  is  very  small.  A  snail  but  not neql ig ible  absorp¬ 
tion  occurs  upon  reflection  at  tie  F2  layer  particularly  near  the  max¬ 
imum  usable  frequency.  Ti.is  absorption  occuiring  at  the  same  time  as 
wave  refraction  isknown  as  devlazivc  absorption.  The  F>  layer  collision 
frequency  is  low  due  to  : iw  molecular  density  so  that  non-devi at i ve 
absorption  within  the  layer  is  nag  I  ig  ible. 

' 4)  For  long  wave  lengths,  tKr  geometric  opt-cs  approximation  is  not 
va  id  particularly  during  daylight.  The  ionosphere  can  be  considered  as 
a  conductor  at  these  frequencies.  A  large  number  of  modes  are  possible 
which  have  the  same  order  of  Intensity. 

b.  Theory.  Equations  for  sk)-wav*»  absorption  have  been  derived  by 
many  theoretical  i  nvest  i  qato  rs .  One  method  of  de-iving  absorption 
eauations,  shown  by  Bremmer  in  " Ts 'rest  rial  Radio  Waves,"  is  indicated 
below. 3  The  geometric— opt ics  approximat ion  fc  r sky-wave  fi  el  d  i ntens i t  ies 
is  appl  •cab1'!  to  an  aosoroiriQ  reu  iv:'  ,  i  i»i  i  1  tjocreit.  In  maosorai'u 

region  the  ind-"x  .if  ref  ract  i'.'i  i'>  ~y—  I  >«  •.  |.  '  ;  •  • • '  ii  i  ft )  form  of  Snell*-* 

Law.  In  tne  natner’atica!  f  i  'r  .1- '  •  n  *■-  -  * >|o-t  'o—Mgiat .  c 
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wave,  i- ex^'ie^t'd!  ter.- appears  accounting  for  tne  atter.uat  ion  per  unit  length 
along  •'(  pat"  t raj ecto'y  .  *ni  s  attenuation  factor  ca"  oe  *r>  t  ten  as- 


Attenuation  Factor 


~  kj  I  (>* )  &  i 

•e  i  * 


Equation 


wh# r«  k  is  the  wave  number,  K  being  the  wave  length 

/R(/i)  is  the  imaginary  part  of  the  refraction  Index 
dl  is  a  differential  distance  along  the  path  trajectory 

When  the  value  for  the  refractior  Index  is  substituted  in  the  equation 
above  and  the  differential  distance  along  the  ray  patn  is  replaced  by 
the  differential  distance  perpendicular  to  the  earth,  this  attenuation 
factor  can  be  written  In  the  following  form:  i  See  3re~mer,  pages  !  76 
a-s  177).  T 

-n  ste  t  r  1  »  . 

- r— 1  J  v(r)  w‘&r 

A  f  -  e  Cu  ri  Equation  2 


where  C  Is  the  velocity  of  light 

r  is  tne  distance  from  the  rent*r  of  the  earth 

n  is  the  number  of  hops 

6  is  the  wave  angular  frequency 

<t  is  the  angle  of  incioence  at  the  ionosphere 

uer  is  the  critical  angular  frequency 

'J(r)  is  the  collisions)  frequency 

The  expression  within  the  integral  sion  Is  independent  of  wave  fre¬ 
quency  for  frequencies  exceeding  the  collision  frequency  This  approx¬ 
imation  is  valid  for  waves  reflected  by  the  £  or  F2  layer  after  passing 
through  the  0  or  E  layer.  The  integral  expression  represents  the  num¬ 
ber  of  collisions  between  molecules  a"d  electrons  situated  in  a  vertical 
colunn  extending  through  the  entire  lower  layer.  it  corresponds  to  non- 
deviative  absorption.  The  expression,  n  sec  <t,  can  be  approximated  by  a 
constant  multiplied  by  the  distance  between  the  receiver  and  the  trans¬ 
mitter.  When  this  is  done  and  .  ' ing  that  all  the  otner  t“rms  in  the 
exponential  expression  can  be  considered  constant,  the  attenuation  fac¬ 
tor  can  be  written  in  Eckersley’s  form: 

A  F  =  e  Equation  3 


where  (J 

f 

(i 


is  a  constant  for  particular  ioniaatior  conditions 
is  the  wave  frequency 

is  the  transmission  distance  between  the  recei.er  and  trans¬ 
mitter 


This  fo  m\  of  absorption  facto'  has  been  used  ir.  the  absorption  calculat- 
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tion  methods  for  long  distance  circuits.  It  should  be  noted  that  the 
effects  of  the  e<irth»s  magnetic  fielc  are  net  included  in  the  formula. 

c.  Afpleton's  Equation.  '  ir  I  ,  J .  Appleton  derived  an  expression 
for  the  total  non—dev  i  at  i  ve  F  la>er  absorption  under  the  condition  that 
the  freqvency  of  the  transmitted  *ave  /'as  n,ur*»  Treater  than  the  colli¬ 
sion  frequency.’*  This  equation  written  sn  terms  of  an  ionosphere  re- 
f,ect:c''  ■''efficient  tor  vertical  inr!*ence  folic v.  s: 


-lot  f.  *  u.lzH  ■ 


‘-■re  z  l  Vn  r  Oa  *• 


Equation  4 


*  '  (*>  ♦  u>  r 


where  f  is  the  reflection  coefficient 
P  is  the  electronic  cfirie 
r  is  the  mess  or  the  elect  rot 
r  is  the  velocity  of  light 

A  q  is  the  electron  density  ’t  the  ■  -i  <  i  n  un  ionization  of  the  layer 
when  the  son's  zenith  .-i  r  -j  I zero 
V0  is  the  co'lihion-l  f-c-u.,.-  .  *t  ’re  height  of  maximum  ioniza- 

t  ion 

H  i s  the  see i e  he i «ht 
V  is  the  zenith  anqle  of  the  sun 

u-  is  the  angular  f  requer  cy  of  tne  :  romga* ed  radio  wave 
& i  is  the  angular  gv ro-magnet ic  frequency  due  to  the  longitudinal 
component  of  the  earf's  Niaqnetic  field 

This  eqjatinn,  based  or.  a  ‘Irpnar  d  ist  i  i  h  .t  ion  ,»f  ionization,  shows  ’tie 
vertical  incidence  aosorpt  i  or  for  the  .-t  .  nary  component  k*  e'ec’to- 
magnetic  wave.  The  abs.  rctitn  /a rip'  no  /.iff.  t  '  \  and  inverse¬ 

ly  with  the  sguPie  of  t1  >.j<  of  t  ’  v  - » •<  'it"  c'..  nd  ’he  I  ong  i  t  urt  i  h,. 1 
comp’onent  of  tne  ay  ro-- 1  i  euuenc  >  . 

d.  Lffect  of  the  E-aptm's  -  i  E  t.»i  upon  ABic’riioN.  Theoretical  eoja- 
tions  show  the  effect  of  the  etmi's  field  .on  absorption.  Tne  rat:o 
of  absorption  in  tne  presence  of  tne  earth's  magnetic  field  to  absorp¬ 
tion  in  absence  of  tne  earth's  naonet  ic  fi-u  is  eqjal  to  tne  ratio  of 
the  imaginary  components  of  '  r-e  res:?"’ive  i  .dices  of  retraction.  iris 
theoretical  ratic,  plotted  1 1  "inure*.  4t>  and  4",  is  given  n>  the  fol¬ 
lowing  equation; 


l;.T<Mi|y=c 


#  f  2  /  —  '  .4 

2  * 


*A  (.-  4  .  :  “*7 


Lq«jat  lor 


where  co  is  the  annular  freouenc, 

Mif  15  the  argular  -iy  *  ;*•  *  r #  :u*'f 
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8  is  the  angle  between  the  earth's  magnetic  field  and  the  radio 
wave 

A  =^.C0B1  8  *-£-8ln*  b 
o>  a 


The  upper  sign  refers  to  the  ordinary  ray  and  the  lower  stgn  to  the 
extraordinary  ray.  This  ratio  for  the  longitudinal  case  (direction  of  the 
earth's  field  parallel  to  wave  normal)  can  be  written  as  foilows: 


When  the  earth's  field  is  perpendicular  to  the  wave  normal,  the  ratios 
for  the  ordinary  and  the  extraordinary  ray  are  given  by  the  following 
expressions: 


•  i  (ordinary  ra;i)  ; 


(extraordinary  ray ) 


e.  absorption  Analysis  introduced.  For  the  absorption  analysis,  the 
effects  on  the  daytime  monthly  median  receiver  input  voltages  of  varia¬ 
tions  in  ray  angle,  running  average  sur.spot  numbers,  gy ro-f requency, 
wave  frequency,  and  the  zenith  angle  of  the  sun  were  studied  in  a  manner 
se'parating  the  individual  variables.  Initially,  the  receiver  input 
voltages  for  each  circuit  and  for  each  month  were  plotted  against  hour 
of  the  day,  as  shown  in  Figure  I.  The  daytime  minimum  of  receiver  input 
voltage  occurs  at  an  hour  corresponding  to  local  noon  at  the  midpoint 
of  the  path  between  the  transmitter  and  receiver  *or  all  circuits. 


f.  Solan  Zenith  Ancle  Relationships.  (I)  fhe  study  of  the  varia¬ 
tion  of  receiver  input  voltages  with  the  sun's  zenith  angle  consisted  of 
determining  the  anqle  at  which  the  daytime  absorption  starts  and  ceases, 
and  the  functional  relationship  between  absorption  cr.c  the  sun's  zenith 
angle.  The ‘sun's  zenith  angle  is  the  anqle  between  the  nerpe-dicu'ar  ti¬ 
the  earth  at  the  ionosphere  reflection  point  and  the  direction  to  the 
sun.  The  cosine  of  this  angle  is  given  by  the  general  formula  of  spher¬ 
ical  trigonometry: 

cos  y  -  stn  £,  sin  L,  *  ccs  L .  cos  I  cos  fl  Equation  6 

where  £(  is  latitude  of  the  ionosphere  reflection  point 
Lt  is  latitude  of  the  sun's  sub-solar  point 
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fl  ig  the  hour  angle  of  the  sun  referred  to  the  ionosphere  re¬ 
flect  ion 

V  is  the  solar  zenith  angle 

In  this  analysis,  solar  zenith  anqles  were  obta ined  graphically  from  con¬ 
tours  of  equal  zenith  angles.  These  contours,  drawn  on  the  same  rec¬ 
tangular  grid  system-used  in  the  Centrai  Radio  Propagation  Laboratory 
"D"  series  predictions  were  computed  by  the  equation  above  with  respect 
to  the  suo-solar  points  for  the  fifteenth  day  of  the  month  .'rom  tne  dec¬ 
lination  of  the  sun  and  the  equation  of  time  given  in  Nautical  almanac 
of  1948.  Since  the  position  of  the  sub-solar  points  do  not  change  much 
from  year-to-year;  the  contours.  Figures  2  through  i3,  can  be  used  for 
any  year. 

(2)  A  study  was  made  of  the  starting  and  stopping  times  of  daytime 
absorption  using  the  curves  showing  monthly  median  receiver  input  volt¬ 
age  plotted  against  xhe  hour  of  the  day.  These  times  occur  when  the 
receiver  input  voltage  curves  begin  to  decrease  at  dawn  and  begin  to 
reach  tfteir  nighttime  value  at  sunset.  The  analysis  for  all  the  cir¬ 
cuits  for  each  month  indicates  that  daytime  absorption  bey  ins  and  ceases 
at  about  the  srme  solar  zenith  angle.  This  angle  is  independent  of 
month,  circuit  length,  and  wave  frequency.  Hgure  M.  shows  the  solar 
zenith  angle  where  absorption  begins  and  where  it  ceases  plotted  against 
wave  frequency.  Figure  16  is  a  mass  plot  showing  thefin?  available  cases 
plotted  over  a  five  year  period.  Figure  17  ts  the  cumulative  distribu¬ 
tion  of  this  mess  plot  data.  The  arithmetic  mean  for  the  starting 
or  stopping  of  daytime  absorption  is'  102.2  degrees.  cigure  Ig  shows  tne 
dawn  and  sunset  time  at  various  heights  above  the  earth.  The  angle  of 
102.2  degrees  indicates  that  daytime  absorption  occurs  between  dawn  and 
sunset  at  an  E  layer  ionosphere  height  of  approximately  115  kilometers. 


(3)  The  functional  relationship  with  the  solar  zenith  angle.in 

Appleton's  absorption  equation  would  oive  zero  absorption  for  an  angle 
of  90  degrees.  Absorption  is  indicated  by  this  analysis  to  scart  or 

cease  at  a  zenith  angle  of  102.2  degrees  and  to  be  maximum  at  a  zenith 

angle  of  zero  degrees.  It  is  o-'ssible  to  set  the  cosine  solar  angle 
relationship  equal  to  zero  at  \.2  degrees  by  multiplying  the  solar 
zenith  angle  by  a  constant  such  that  the  product  of  the  constant  and 
102.2  would  give  a  value  of  90  degrees.  Thus,  absorption  would  start 
or  cease  at  tne  observed  angle  and  would  be  maximum  when  the  sun  Is 
overhead.  Based  on  the  above,  the  functional  relationship  of  absorption 
with  the  angle  was  assumed  to  be  of  the  following  form  in  this  report: 

(COS  app-  (cos  . 881 *PP  Equation  7 

where  a  is  a  constant  equai  to  90°/l02 .2°  =  •■881 

»  is  •  constant  to  be  obtained  from  the  data 

*p  is  the  solar  zenith  angle 

(4)  The  monthly  median  receiver  Input  voltages  were  plotted  for  all 
circuits  for  each  month  against  the  cosine  of  .881  as  shown  In  Figure 
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18.  In  many  cases  the  morning  curve  is  not  as  symmetrica!  as  tne  aftar- 
noon  cy.ve  due  to  lack  of  ionospheric  support  during  the  morning  hours. 
A  composite  curve,  Figure  19,  is  drawn  between  the  valies  for  morning 
and  the  values  for  afternoon.  From  this  composite  curve,the  voltage 
value  at  COS  .881  equal  to  zero  was  taken  as  a  reference.  In  cases 
where  this  value  was  not  obtainable  from  the  data,  the  reference  volt¬ 
ages  were  obtained  hy  extrapolating  the  curves  to  COS  .681  Y-  equal  to 
zero.  Absorption  values  were  obtained  as  the  decibel  difference  between 
the  reference  and  voltage  values  at  particular  solar  angles.  These  ab¬ 
sorption  values  were  then  plotted  against  COS  .881  on  log  paper  as 

shown  in  Figure  20.  These  curves  are  nearly  straight  lines.  Tne  ab¬ 
sorption  curves  described  above  are  not  exactly  useful  since  absorption 
is  also  a  function  of  running  average  sunspot  number  and  the  angle  of 
incidence  at  the  ionosphere.  In  order  to  separate  these  variables,  it  is 
desirable  to  convert  oblique  absorption  values  to  vert ical  incidence  at  a 
particular  sunspot  number.  The  following  paragraphs  indicate  how  the 
oblique  incident  absoiption  values  were  converted  to  vertical  incidence 
values  for  sunspot  number  zero. 

g.  Solar  Activity  Effects.  In  order  to  study  the  variation  of  day¬ 

time  ahsorptlon  with  average  sunspot  number,  it  was  necessary  to  analyze 
data  coverlnq  a  period  of  several  years.  cigure  21  shows  absorption  ard 
12  month  running  average  sunspot  number  plotted  against  time  for  four 
routes.  There  was  an  initial  attempt  to  analyze  the  absorption  for  a 
particul ar  value  of  the  sun's  zenith  angle.  A  higher  correlation  coef* 
ficient  was  obtained  for  the  analysis  of  the  absorption  at  midday. 
Twelve  month  running  average  mi.dday  absorption  was  plotted  against  12 

month  running  average  sunspot  number  for  the  eight  routes  tnat  had  data 

over  a  sufficiently  long  time  period.  The  circuits  analyzed  cover  vari¬ 
ous  transmission  distances  from  pa  kilometers  to  3636  kilometers  and 
frequencies  from  about  2  megacycles  to  10  megacycles.  By  plotting  12 
month  running  averages,  the  seasonal  trends  of  absorption  were  removed. 
The  data  shows  that  absorption  varies  linearly  with  sunspot  number  as 
shown  in  the  following  equation: 

-  And  +l,a!  equation  8 

where  3  is  the  runninq  average  sunspot  number 

b  is  a  constant  obtained  from  the  experimental  data  equal  to  .0037 
A0  is  the  absorption  at  Sunspot  Number  o 

Az  is  the  absorption  at  Sunspot  Number  X 

The  slope  of  the  linear  relation  between  absorption  and  sunspot  number 
is  .0037  as  an  average  for  the  eight  circuits.  Table  2  lists  the  values 
of  the  slope  o.d  the  linear  correlation  coefficient  obtained  from  the 
least  sevares  ana'ys's. 

h.  LAYER  HEIGHT  VARIATIONS.  t|!  F.  VY.  G.  White  mentions  that  the 
variation  of  absorption  between  winter  and  summer  i  s  less  than  that  expect¬ 
ed  by  theory  based  onthesolar  zenith  angle.4  Figure  22  shows  the  midday 


3 


absorption  ext  *apo !  ated  *0  a  aero  Jeoiee  angle  for  the  radio  cir¬ 

cuit  between  f/exico  City  and  Boston.  *'..s  f  i  ou  re  shows  the  minimum  F2 
layer  virtual  heights  observed  -t  c«  u  !  "  1  ana  at  ate  'Jrivprsity  for  me 
same  period  uf  time  covered  by  th“  absorption  -".jrves.  Tms  minimum  vir¬ 
tual  height  curve  is  plotted  with  intreisirq  height  in  a  downward  direc¬ 
tion.  By  com.parinq  th  •  <ro*>th— to-mort*.  itr  at  ions  of  midday  absorption 
with  the  monthly  average  co  'ayei  w'd-lav  heig.t,  it  .s  seer  that 

the  absorptior.  appears  to  losio.*.  1  ?  '<\vf-r  heights  nverse.y.  Curing  the 
summer  season,  the  rp  '  ;vzr  r.e'ghts  am  iqn-r  and,  .hus,  the  ray  angles 
at  the  ionosphere  are  steeper  tt.m  <>ie  w  ater  scasrr,.  Ihi. 

tends  to  decrease  absorptior.  In'irg  tie  sw'et  -oaths  relative  the 
values  during  the  winter  months,  explaining  the  difference  between  uin- 
ler  and  summer  absorption  not  accountable  c\  tne  zer  th  angle  alone. 

(2)  The  diurnal  variation  of  minimum  monthly  median  F?  layer  virtua1 

heights  to'  earh  mont*-  are  shown  in  r.yures  21  tl  rough  These  fio- 

ures  reoresent  the  average  F2  liycr  vi'tual  teigv.  at  all  ionosphere 
stations  as  reported  by  the  batioral  Bureau  cf  otancards  in  ihe  Central 
Nadio  Propagation  Laboratory  f  Leriot  '-ports  for  'O45  through  1943. 
Based  on  this  analysis,  f2  layer  v'rtual  he  yhts  do  not  appear  to  have  a 
regular  variation  pattern  throughout  the  ->o'ar  cycle  or  between  geo¬ 
magnetic  zones.  Figures  2!-  through  J4  may  therpCore  be  used  to  approxi¬ 
mate  the  minimum  F2  iayer  virtual  neiaht  it  nr.y  geographic  location  at 
any  time.  Figure  3£  shows  that,  in  so-n  ! ©.-at  :ons  the  mini'um  ??  'ayer 
virtual  height  increases  whi  I  e  in  others  ,t  fpfrpisi".  or  lemams  constant 
as  the  moving  average  sunspot  nurne'  me  rtas.-s,  rigu'o  Jf>  shows  the 
average  heights  obtained  f  rom  the  contours  and  the  observed  neon  f2  layer 
heiahts  at  Washington,  C.  C.  for  five  years.  ’he  monthly  median  F  laver 
virtual  heights  do  not  vary  oreatlr  from  month-to-morth.  iicorj  }~  shows 
the  minimum  virtual  £  layer  heights  ar  observed  in  A»>'“  -ct.n,  L\  c.  i  n 
June  and  December  of  194*  plotted  ac.'.irq  vi.ir  of  tie  day. 

(3)  Only  F2  layer  a. id  £  layer  transmission  modes  are  considered  1 " 
this  report.  E  layer  transmission  mono-*  a"  acs;-ed  to  be  'ef'erfed 
from  a  110  kilometer  virtua'  heisnt.  1 ;n  r  \  ransm.i  r  s  i  modes  *cr 
equai  antenna  radiations  probably  I'd  between  ihe  Fy  and  tne  fc  mod’s  in 
intensity.  General'y,  rays  that  penetrate  the  1  layer  w..'  be  reflected 
by  the  F2  layer  and  Ft  iayer  mode  iiubat.1.  are  s~a.  1-'  ’  w  rags  tude 
than  the  F2  layer  modes. 

1.  Conversion  to  vertical  ncimnlc  ’  hbr.orrt.cn  was  assumed  to 

be  directly  proportional  to  tne  ar.s.e  •  ecia<.ncr  at  a  .PC  k.iomeiei 
height.  These  angles  were  based  ,n  oir  'e<'ecte'1  ra.»s  between  ? 
curved  earth  ar.d  a  curved  ionospne'e.  ;  •  t  venieice  in  analysis,  tne 
secant  ot  the  angle  at  ■P''1  r  .  tcme*..-.  r.»  aten  to  p?  anni«>  between 

the  ray  and  ground.  figure  -?  si-.,*  -n<  - .... t  of  the  art:  e  of  inci¬ 

dence  for  the  various  'aye'  height'-  *s.*  ‘  r  •'  i*.  pi  g'ej-  ci"-<.  distj"ce 
and  angle  between  the  ray  a-.i  the  g  -.  1  »■  -cca-'t  factors  f-'r  muiti- 

hop  modes  are  eaua!  to  the  seemt  fa-.*  *  *0*  s.rq  •'  hup  r.-'  !(s  of  the 
same  radiat:on  angle  n't  p'.eq  ■  •  .*  '  ‘  on-sphe'e  *ef  lec  <  ’ 
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For  tne  various  possible  propagation  modes  reflected  from  the  F.  ^nd  F 2 
layers,  the  secant  factors  at  a  100  kilometer  heigh?  were  obtained  at 
the  radiation  angles  determined  by  the  number  of  ionospheric  reflections, 
the  vi.t-ial  layer  height,  and  transmission  distance.  Maximum  usable 
frequency  curves  were  plotted  for  the  transmission  path  for  the  period 
of  the  data  from  the  observed  critical  frequencies  shown  in  the  Central 
Radio  Propagation  Laboratory  F  Series.  These  curves,  see  Figure  39,  show 
which  sky-wave  modes  were  possible  at  certain  hours  and  frequencies. 

(2)  The  oblique  incidence  absorptions  were  converted  into  vertical 
incidence  absorptions  by  dividing  by  the  serant  factors.  The  layer 
heights  corresponding  to  the  particular  loca!  time  at  the  ionosphere 
reflection  point  were  determined  from  Figure  23  through  34,  which  give 
the  average  layer  heights  for  each  month  of  the  year.  Since  both  the 
solar  zenith  angles  and  thj  layer  heights  were  obtained  frcm  charts  in 
terms  of  latitude  and  local  time,  the  secant  fa^lo's  for  F2  layer  modes 
were  obtained  at  corresponding  solar  2e. ith  angles.  In  addition  to 
dividing  the  oblique  incidence  absorptions  by  the  secant  factors,  they 
were  divided  by  the  solar  activity  factor  given  in  Equation  8  to  convert 
the  vertical  incidence  values  to  sunspot  number  zero. 

(3)  The  vertical  incidence  absorption  at  sunspot  number  zero  for  each 
momh  and  for  ail  the  c!rr.uits  was  then  plotted  against  the  cosine  of 
.68'^"  Figure  40  shows  a  typical  plot  of  this  kind  end  indicates  the 
mathematics  used  it  obtaining  the  vertical  incidence  absorption  for  sun- 
soot  number  zero  from  the  oblique  incidence  absorption  at  a  Known  sun¬ 
spot  number.  These  curves  are  more  linear  than  the  plots  of  oblique 
incidence  absorption  against  the  same  function  of  solar  zenith  angle. 
Secant  factor  variations  due  to  ionospheric  layer  height  variations  are 
accounted  for.  in  this  analysis,  no  account  is  made  of  the  chai  -e  !n 
antenna  r -d i at i on  due  to  hourly  changes  in  iayer  heights.  It  is  fe.  • 
that  apparent  absorption  variations  due  to  differences  in  antenr.a 
radiations  that  occur  •from  hour-to-hour  do  not  cause  large  errors  when 
all  the  routes  arc  considered  as  a  group.  Some  of  tne  circuits  have 
higher  total  antenna  gains  at  lower  angles  than  at  higher  angles,  while 
other  ci rcuits  have  opposite  characteristics. 

(4)  The  plots  of  vertical  incidence  absorption  against  the  cosine  of 
.88I<A  were  extrapolated  to  the  value  of  cosine  .88I<1'  equal  to  I,  the 
sub-solar  absorption  forvertical  inc  idence  at  sunspot  number  zero.  These 
values  should  not  vary  with  season  since  the  effect  of  height  variations 
have  been  accounted  for.  Table  3  shows  that  residual  seasonal  fluctua¬ 
tions  ere  negligible.  The  effect  of  transmission  distance  upon  these  ver¬ 
tically  incident  values  should  also  have  been  removeo.  Figure  4i  shows 
the  verticel  incidence  absorption  plotted  against  transmission  distance 
for  those  values  of  frequency  for  which  data  is  available  at  several 
distances.  No  appreciable  variations  were  noted  with  transmission  ois- 


tance.  The  exponent  of  the  relationship  between  vertical  incidence 
absorption  and  COS  .881  4  varied  somewhat  between  circuits.  These 
exponents  are  plotted  as  a  mass  plot  against  the  months  of  the  year  in 
Figure  42.  Figure  43  shows  the  cumulative  distribution  of  the  exponents 
for  all  the  observations  from  the  data,  plotted  on  normal  probability 
paper.  Figure  44  shows  the  exponents  plotted  against  wave  frequency. 
Table  3  Is  a  summary  of  the  seasonal  variations  of  the  exponents.  It  is 
noted  that  the  seasonal  average  values  are  very  nearly  alike  and  that  the 
composite  median  for  all  available  observations  is  1.30,  (Figure  431. 

15)  For  distances  greater  than  4000  kilometers,  the  conversion  of 
oblique  incidence  absorption  to  vertical  incidence  was  performed  in  a 
slightly  different  manner.  The  averaqe  of  the  layer  heiahts  between 
2000  kilometer  points  from  the  transmitting  and  receiving  locations  was 
used  in  determining  the  radiation  angles  above  the  ground  plane  for 
various  probable  transmission  modes.  For  these  modes,  the  solar  zenith 
angles  were  determined  at  points  corresponding  to  the  positions  of  iono¬ 
sphere  reflections.  The  average  of  the  cos  .881  't-  for  these  reflection 
points  was  taken  as  the  value  for  the  circuit.  The  oblique  incidence 
absorptions  were  divided  b j  secant  factors  for  the  most  probable  mudes 
and  by  the  solar  activity  facto'-  to  transform  the  values  to  vertical 
incidence  and  sunspot  number  zero.  Figure  45  shows  the  secant  factors 
at  a  too  kilometer  height  plotted  against  distance  for  various  number 
of  hops  for  a  reflecting  layer  height  of  300  kilometers.  Figure  38  was 
used  to  obtain  the  secant  factors  in  this  report,  the  secant  factors  for 
multi-hop  modes  being  equal  to  that  for  single  hop  modes  of  the  same 
radiation  angle  multiplied  by  the  number  of  hops.  The  vertical  incidence 
absorption  values  at  sunspot  number  zero  were  then  clotted  against  the 
average  cosine  of  .881  V  as  in  the  case  of  shorter  routes.  The  subsolar 
absorption  values  at  sunspot  verc.  and  s lopes  are  included  with  the  values 
for  shorter  routes  in  Figures  41,  42.  43,  and  44. 

j.  OasEkvto  Effects  of  'he  Eanth-s  Magnetic  Field,  iii  Since  ab¬ 
sorption  is  a  function  of  both  wave  frequency  and  gyro-frequency,  it  is 
desirable  to  analyze  the  relationship  between  absorption  and  wave  fre¬ 
quency  in  a  manner  eliminating  the  effects  of  the  magnetic  fieid.  The 
Appleton  equation  for  the  ordinary  ray  indicates  that  absorption  varies 
inversely  with  the  sauare  of  t  he  sum  of  the  wave  frequency  and  the  lono- 
itudinal  component  of  the  gy ro- Frequency.  theoretical  equation  5  snows 
iiun  absorption  values  can  be  obtained  in  absence  o*  the  earth's  magnetic 
field.  These  absorption  values  should  show  a  consistent  relationship 
with  wave  frequency  since  the  magnetic  field  effect  has  been  removed. 
In  this  report,  the  relationship  bstween  absorpt'on  and  wave  fre.uency 
was  studied  in  absence  of  the  magnetic  field  and  by  a  more  approximate 
method. 

(2)  Figures  46  and  47  can  be  used  to  --op.vert  absorption  to  values  in 
absence  of  the  earth's  field.  These  curves  show  the  relatioiship 
between  absorption  with  and  without  the  presence  of  the  earth's  field 
for  both  the  ordinary  and  the  axtraord inary  rays,  in  order  to  use  these 


curves,  the  angle  between  the  earth's  magnetic  fiejd  and  the  ray  direc¬ 
tion  must  te  known  in  addition  to  the  gy ro-f requency . 

(3)  Two  space  angles  in  a  spherical  coordinate  system  determine  the 
direction  of  any  line.  The  variation  and  tne  dip  determine  the  direc¬ 
tion  of  the  earth's  field.  Figure  48  snows  the  variation  of  the  earth's 
field  while  Figure  49  shows  its  dip.  The  angle  of  incidence  at  a  IOC 
kilometer  height  and  the  bearing  of  the  great  circle  path  between  the 
transmitter  and  receiver  at  the  ionosphere  reflection  point  determine 
the  direction  of  a  ray  at  the  absorption  height.  The  variation  and  dip 
of  the  earth's  field  and  the  bearing  of  the  great  circle  path  are  deter¬ 
mined  graphically  by  using  an  overlay  and  charts  drawn  on  the  National 
Bureau  of  Standards  rectangular  grid  system.  Figure  bO  is  a  bearing 
chart  which  has  been  devised  at  the  Radio  Propagation  Unit  to  determine 
tne  bearing  from  north  for  any  point  along  the  great  circle  path  between 
any  two  points  on  tne  earth's  surface. 

(4)  a  transparent  piece  of  paper  is  placed  on  a  world  map,  Figure  83, 
and  the  positions  of  the  transmitting  and  receiving  stations  are  marked 
in  addition  to  drawing  a  line  over  the  equator.  This  overlay  is  then 
placed  on  tne  great  circle  chart.  Figure  84,  and  moved  horizontally  with 
the  equatorial  lines  coinciding  until  the  transmitting  and  receiving 
locations  tie  on  the  same  great  circle  or  proportional  distance  between 
adjacent  great  circles.  The  great  circle  is  sketched  on  the  transpar¬ 
ency  to  the  point  where  the  great  circle  path  crosses  the  equator.  The 
transparency  is  tnen  placed  on  the  Radio  Propagation  unit  Bear;";  Chart 
with  the  equatorial  lines  coinciding  and  the  point  where  the  great 
circle  pat-n  crosses  the  equator  is  placed  on  one  of  the  two  equatorial 
polar  points.  The  bearings  are  found  in  either  direction  at  any  point 
along  the  q.-eat  circle  path.  It  a  point  on  the  great  circle  lies  over  a 
section  of  the  bearing  chart  marked  E  (or  W)  the  bearing  is  found  for  an 
easterly  (westerly)  direction.  The  overlap  is  then  placed  over  the  dip 
and  variation  charts  to  find  these  values  at  the  path  midpoint.  The 
over'ey  is  also  used  to  find  the  gy ro-f requency  at  the  ionosphere  re¬ 
flection  point.  Figure  ol  shows  tne  gy  ro-f  requency  chart  drawn  on  the 
same  rectangular  qrid.  The  values  are  obtained  from  the  measured  total 
earth's  magnetic  field  at  the  surface  from  the  Hydrographic  Office 
Charts  O'*  194b  extrapolated  to  a  100  kilometer  height  by  Schmidt's 
equat ion.5 

(b)  Since  the  directions  of  the  earth's  field  and  the  ray  path  are 
known,  the  anqle  between  the  two  can  be  determined  by  the  general  for¬ 
mula  of  spherical  trigonometry  or  from  the  D'Ocagne  Nomogram,  Fiqure  ai. 
All  the  variables  reauired  for  the  determination  of  *iie  factor  resating 
absorption  in  presence  of  tne  earth's  field  with  absorption  in  the  ab¬ 
sence  of  the  earth's  field  are  known  since  the  angle  between  the  ray 
direction  and  the  earth's  field  is  known  in  audition  to  the  gyro-frequency 
and  the  wave  frequency.  In  this  analysis,  separate  factors  are  deter¬ 
mined  for  both  entry  and  exit  since  the  angle  between  the  earth's  field 
and  ray  differs  for  the  upooing  and  downceming  parts  ot  the  ray  paths. 
The  average  of  the  factors  for  entry  and  exit  were  used  to  corverl 
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absorption  to  that  expected  .vhcn  the  earth's  field  is  not  present. 
Figure  53  shows  this  absorption  in  absence  of  the  earth's  magnetic 
field  at  vertical  incidence  plotted  against  the  wave  frequency  tor 
ordinary  ray  transmission.  The  curve  includes  absorption  values  tor  all 
four  seasons.  Monthly  mocian  values  for  each  individual  circuit  are 
shown  by  dots. 

(6)  Due  to  the  complexity  required  in  determining  the  angle  between 
the  earth's  field  and  the  geometric  rays,  it  is  desirable  to  use  s 
method  which  closely  approximates  the  effect  of  the  earth's  field  but 
does  not  require  the  determi nat ion  cf  the  angle.  Various  approximations 
were  i rvest (gated.  The  data  indicates  that  wnen  absorption  is  plotted 
against  the  sum  of  the  wave  frequency  and  the  total  gyro-frequency  that 
a  consistent  relationship  is  obtained.  F’gure  34  shows  the  equivalent 
vertical  absorption  for  sunspot  number  zero  and  solar  zenith  angle  zero 
plotted  against  the  sum  of  the  wave  frequency  and  the  gy ro-f requency  for 
all  four  seasons  of  the  year.  The  individual  dots  represent  monthly 
median  values  fer  each  circuit.  Figure  53  shows  the  same  curve  when 
the  yearly  average  values  are  used.  Least  squares  straight  lines  based 
on  the  yearly  average  values  are  drawn  on  a|l  these  plots.  The  slope 
of  these  lines  shows  that  absorption  varies  closely  in  an  inverse 
squared  relationship  with  the  sum  of  the  wave  frequency  and  the  total 
gyro-frequency.  It  is  noted  that  few  experimental  points  were  available 
below  about  5  megacyles  and  ,  thus,  the  reliability  of  the  curve  is  less 
at  these  lower  frequencies.  Thus,  the  data  shows  that  vertical  incidence 
absorption  for  a  solar  zenith  angle  of  zero  degrees  and  a  sunspot  number 
of  zero  can  be  related  to  the  sum  of  the  wave  frequency  and  the  gyro- 
f requency  by  the  following  equation: 


Frequency  Function  = 


Corn  tan t 

(f  +  f/ 


Equation  9 


where  ft  is  the  gyro-frequ  nCy 

f  is  the  wave  frequency 
d  is  a  cor-.tant 

k.  Radio  Propagation  unit  Equation.  The  statistical  analysis  of  the 
monthly  median  receiver  input  voltages  indicates  that  absorption  fits  a 
formula  similar  to  Appleton's  converted  to  obiique  incidence. 


A 


in  decibels  =  2L*2£JL  HUL 


Equation  10 


where  is  the  sun's  zenith  angle 
Q  is  a  constant  =  .881 
b  is  a  constant  =  .0087 
C  is  a  constant  -  61b. 0 
«  is  a  constant  =  1.8 
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'i  is  a  constant  = 1.98 
f  Is  the  wave  frequency 

f3  is  the  qy ro-f requency  at  a  IOC  kilometer  height 
n  is  the  rumber  of  ionospheric  reflections 

<p  is  the  angle  at  a  100  kilometer  height  between  a  perpendicular 
to  the  earth  and  the  ray  path 
S  is  the  12  month  running  average  sunspot  number 

It  should  be  noted  in  the  above  equation  that  the  total  gyro-frequency 
is  used  instead  or  the  longitudinal  component  and  that  this  gyro- 
frequency  is  the  value  at  a  100  kilometer  height. 

3.  Analysis  or  Statistical  Variation  of  Hourly  Median  Sky-wave  Field 

1  ntens i ties. 

a.  General.  The  mass  plots  of  the  hourly  median  receiver  input 
voltages  were  used  in  studying  the  day-to-day  variations  over  monthly 
and  yearly  time  periods.  Figuru  64-  shows  a  typical  mass  plot.  The 
dots  show  the  number  of  days  that  the  hourly  median  voltages  lie  within 

2  decibel  class  intervals  over  the  month.  Separate  studies  were  per¬ 
formed  for  day  and  night  hours.  These  studies  indicated  that  when  a 
sufficiently  long  tine  period  was  used  both  the  nighttime  and  the  day¬ 
time  day-to-day  variations  were  similar  for  each  circuit  and  are  prac¬ 
tically  independent  of  wave  frequency,  transmission  distance,  ana  nour 
of  the  day. 

b.  Nighttime  Day-to-day  variations.  In  the  analysis  of  nighttime 
voltages,  the  statistical  study  was  made  at  the  night  hour  at  which  the 
monthly  median  value  was  the  highest.  The  data  was  plotted  in  histo¬ 
grams,  indicating  the  number  of  cases  wnere  the  voltage  fell  between  4 
decibel  intervals,  as  shown  in  Figure  63,  Cumulative  distribution  plots 
were  made  for  each  month  for  all  the  circuits.  These  distributions  were 
plotted  both  on  normal  probability  paper  and  on  Rayleigh  coordinate 
paper,  as  snown  .1  Figures  66  ana  67.  The  number  of  cases  within  4 
decibel  class  intervals  were  also  plotted  against  too  months  of  the 
year.  Figure  68  shows  one  of  these  plots  »'hich  were  used  in  analyzing 
the  yearly  trend  of  t..e  monthly  median  values.  For  each  circuit,  the 
cumulative  sums  of  ati  the  enses  within  4  decibel  class  intervals  for 
the  entire  year  were  plotted  on  Rayleigh  coordinate  graph  paper.  These 
curves  are  shown  for  each  route  of  this  report,  Figuies  69  through  72. 
Comparing  the  v»rious  figures,  it  can  be  noted  that  they  are  nearly 
linear  between  percentages  exceede.  90  percent  of  time  and  10  percent 
of  time.  Receiver  input  voltage  in  decibels  referred  to  the  yearly 
median  were  'plotted  against  percentage  of  the  time  the  voltage  is  ex¬ 
ceeded.  Figure  73  is  a  cumulative  distribution  that  includes  all  the 
circuits.  The  decibel  values  referred  to  the  median  for  each  circuit 
were  read  from  the  cumulative  distribution  plots  at  particular  percent¬ 
age  values  and  these  W6r6  plotted  as  a  composite.  The  composite  is 
nearly  a  straight  line  with  a  slope  approximating  that  obtained  by  plot¬ 
ting  a  Rayleigh  distribution  on  this  paper.  There  are  approx imate I y  13 
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decibels  between  the  exceeded  10  oercent  of  the  time  and  the  exceeded 
90  percent  of  the  time  values.  Figure  74  shows  the  cumulative  distri¬ 
bution  of  the  decibel  differences  between  the  exceeded  10  percent  of  the 
time  and  the  exceeded  90  percent  of  the  time  values  for  each  circuit 
when  tnese  are  taken  on  a  yearly  basis.  This  figure  is  useful  in  show¬ 
ing  the  percentage  of  circuits  having  particular  decibel  differences 
between  the  10  percent  exceeded  and  the  90  percent  exceeded  values. 

r.  Daytime  Day-to-day  variations.  (I)  The  daytime  voltages  were 
analyzed  in  a  similar  manner.  Ihese  again  show  considerable  differences 
if  individual  routes  are  examined  on  a  monthly  basis.  But,  again,  if 
the  data  for  a  pariod  as  long  as  a  year  is  used  in  analyzing  the  day-to- 
day  variations,  the  routes  appear  simi'a-.  For  the  daytime  /oltages, 
there  were  two  separate  analyses,  in  one,  the  data  was  analyzed  at  four 
six  hour  intervals;  6  a.  m.,  noon,  6  p.  m. ,  and  midnight.  For  each 
circul*  at  these  nours,  the  decibel  differences  between  the  monthly 
median  and  percentages  exceeded  10,  30,  i 0,  and  90  percent  of  tha  time 
were  obtained  from  the  mass  plots  by  counting  the  number  of  cases  above 
and  below  the  monthly  median  values.  The  average  differences  between 
these  Dercentaaes  for  all  the  circuits  at  tne  four  hours  of  the  day 
were  combined  into  composite  plots  at  approximately  the  same  wave  fre¬ 
quencies.  Figure  7b  shoes  this  composite  cumulative  distribution  for 
particular  wave  frequency  values. 

(II  'n  .he  other  analysis,  the  variations  wa re  obtained  at  a  sol*r 
zenith  angle  of  60  degrees.  The  zenith  anqle  of  60  degrees  was  chosen 
to  have  data  throughout  the  year  for  a  large  number  of  circuits.  Large 
zenith  angles  are  not  possible  at  al I  locations  for  the  entire  year.  The 
receiver  input  voltages  were  grouped  into  2  occibel  class  intervals  for 
for  both  morning  and  afternoon  fer  each  month  for  all  the  circuits. 
Tigure  76  is  a  typical  histogram  snowinq  the  number  of  instances  at  2 
decibel  class  intervals  referred  to  the  monthly  median  for  both  morning 
and  afternoor  at  a  solar  zenith  angle  of  60  degrees  for  circuit  for 
one  month.  The  cumulative  distributions  were  calculated  and  plotted. 
From  the  plotted  cumulative  distribution  curves  for  each  circuit  and 
each  month  for  morniig  and  afternoon,  the  decibel  differences  between 
the  median  value  were  read  at  percentages  exetode:  !,  b,  10,  3C,  70, 
90,  9b,  and  99  percent  of  the  time.  These  .-e.c  averaged  for  all  the 
months  of  tne  year  w:th  betn  morning  and  afternoon  values  included. 
Figures  77  and  78  are  cumulative  distributions  for  each  circuit  on  this 
yearly  basis.  Figure  79  is  a  composite  for  all  the  circuits.  The  com¬ 
posite  curvp  for  all  the  circuits  during  daylight  hours  again  plo*s 
nearly  linearly  on  Rayleigh  coordinate  paper.  The  aecibe!  difference 
between  tho  exceeded  10  and  90  percent  cf  the  time  values  is  14  deci¬ 
bels,  nearly  that  expected  for  a  Rayleigh  d i st r i but  ion .  The  decibel 
differences  between  the  10  and  90  percent  values  were  plotted  for  all 
the  routes  'or  the  solar  zenith  ancle  for  60  degiees  as  a  cumulative 
distribution  on  Figure  80. 

d.  Summary.  (I)  Smith,  and  harnington  in  Central  Radio  Propagation 
Laboratory  Report  i  -  6  show  that  the  distribution  of  day-to-day  hourly 
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median  field  intensities  over  long  time  periods  is  roughly  ncrmqi  in 
decibel*,  between  exceeded  10  to  90  percent  of  the  time  value*..  They 
show  a  12  decibel  difference  between  the  exceeded  10  anu  the  exceeded 
90  percent  of  the  time  values  as  a  median  for  the  circuits  analyzed. 
The  standard  deviation  of  the  approximate  normal  distribution  was  shown 
to  be  about.  t>  decibels.  Figure  61  shows  a  normal  distribution  with  a 
6  decibel  standard  deviation  plotted  on  Rayleigh  coordinate  paper. 
Short  term  variations  in  sky-wave  fie'd  intensities,  appear  to  follow 
tne  Rayleigh  distribution,  given  by  tne  rol lowing  equat'on*. 

P  =  100e~'"31fc)  Equation  || 

where  E  is  the  instantaneous  field  intensity 
Ew  is  the  median  field  intensity 

P  is  the  percentage  of  the  time  the  instantaneous  value  is  exceeded 

(2)  as  indicated  in  the  preceding  paragraphs,  both  daytime  and  night¬ 
time  hourly  median  receiver  input  voltages  plot  nearly  linearly  on 
Rayleigh  coordinate  paper  with  abouf  the  same  slope  as  expscted  from  a 
Rayleigh  distribution.  Thus,  the  day-to-day  variations  over  a  period 
as  long  as  a  year  have  approximate  Rayleigh  distributions.  A  Rayleigh 
distribution  is  the  same  as  the  distribution  for  the  magnitudes  of  a 
complex  variate  with  equal  standard  deviations  along  two  rectangular 
axes.  The  work  of  Ray  S.  Hoyt  of  the  Bell  System  on  probability  theory 
for  two  dimensional  normal  distributions  can  be  cited.*  In  his  i 933 
paper,  he  shows  cumulative  distributions  for  the  magnitudes  for  various 
values  of  a  function,  b. 

1  r  (sv/  su) 

b  «  - - - Equation  1 2 

'  +Msu) 

where  Sp  is  the  standard  deviation  along  the  U  axis 

Sj  is  the  standaro  deviation  along  the  V  axis 

When  the  standard  deviations  along  two*  rectangular  axes  are  equal,  the 
value  of  the  function  is  equal  to  zero.  When  one  standard  deviation 
greatly  exceeds  the  other,  the  function  is  equal  to  I.  Cumulative 
distributions  for  b  values  of  0  and  I  are  also  shown  on  Figure  81. 

(31  The  day-to-day  variations  over  as  short  a  period  as  a  month 
vary  from  month-to-month  and  between  circuits,  Many  of  these  variations 
can  be  accounted  for  by  variations  in  ionospheric  support.  Figures 
showing  the  critical  frequencies  for  transmission  by  various  sky-wave 
modes  are  usaoie  in  analyzing  reasons  for  many  of  these  variations. 
The  extent  of  ionization  of  the  layers  of  the  ionosphere  oetermines 

which  modes  are  possible.  Many  of  tne  variations  can  be  traced  to  solar 
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activity.  Figure  82  shows  the  nighttime  monthly  median  receiver  input 
voltages  for  MiT  receiving  XEW.V  at  two  night  hours.  Night  values  remain 
high  between  early  1946  and  midyear  1948,  the  peak  of  the  present  solar 
cycle.  This  is  probably  due  to  the  fact  that  the  ionosphere  is  suffi¬ 
ciently  ionized  to  consistently  support  the  F2  two  hop  transmission 
mode,  the  predominate  sky-wave  mode  for  that  circuit. 

4.  The  Analysis  of  the  Absolute  Magnitude  of  Field  Intensities 

a.  General.  The  ■■.bsolute  magnitudes  of  the  sky-wave  field  intensi¬ 
ties  were  analyzed  ov  using  the  measured  data  and  circuit  information. 
Theoretical  receiver  input  voltages  were  compared  with  observed  data 
values  in  microvolts  receiver  input.  Measured  instantaneous  output 
voltages,  recorded  on  cont  i  nuous  tste  rl  i  ne-Angus  type  recorders,  were 
calibrated  in*terms  of  receiver  input  voltage  at  the  receiving  stations. 
Hourly  median  values  were  computed  at  the  rcceivinq  locations  *rom  the 
instantaneous  recordings.  I nformat ion  about  the  t ransmi tt ing  and  receiv¬ 
ing  antenna  gains,  the  transmitting  antenna  radiated  powers,  and  the 
receiving  circuit  input  characteristics  was  obtained  directly  from  the 
various  transmitting  and  receiving  stations.  Measured  antenna  gains 
were  available  for  a  few  of  the  cirruits.  The  antenna  gains  were  calcu¬ 
lated  in  this  analysis  from  standard  equ  tions  such  as  shown  in  Radio 
Propagation  uni't  Report  No.  7  for  linear  thin  wire  antennas,  accounting 
for  imperfect  earth  conditions  ano  assumed  antenna  losses.®  Table  4 
shows  the  transmitting  antenna  gains  and  the  receiving  antenna  effective 
heights  for  various  reciat:on  angles  above  the  ground  plane.  Mismatch 
losses  based  on  ihe  input  circuit  characteristics  were  included  in  the 
receiving  antenna  effective  heights. 

b.  the  Calculated  Receiver  Input  Voltages.  Mil  The  measured  and 
calculated  receiver  input  voltaqes  were  compared  for  niaht  hours.  For 
daytime  hours,  the  absolute  magnitude  can  be  obtained  from,  the  niqht 
values  by  accounting  frr  the  daytime  absorption  losses.  The  receiver 
Input  voltages  were  calculated  for  various  possible  ray  paths  for  assumed 
lossless  reflections  at  the  ionosphere.  These  values  show  what  the 
receiver  input  voltage  would  be  if  the  only  losses  are  those  due  to 
free  space  distance  attenuation  and  ground  reflections. 

(2)  The  free  space  loss  along  a  ray  path  is  dependent  upon  distance 
between  the  transmitter  and  receiver,  the  v'rtual  layer  height,  and  the 
number  of  hops.  These  ware  calculated  for  various  assumed  layer  heights. 
The  analysis  indicates  that  within  a  very  close  approximation  that  the 
ratio  between  the  distance  along  a  ray  path  and  the  distance  aiong  the 
curved  earth  between  the  transmitter  and  receiver  is  indeoendent  of 
layer  height  for  equal  radiation  angles.  Thus,  tne  tree  space  field 
intensity  ot  any  ray  can  be  given  ns  a  function  of  tne  radiation  onqle 
and  the  transmission  distance  alo^e.  The  incident  field  intensity  in 
free  space  fora  reference  inverse  distance  field  intensity  can  he  qiven 
by  the  following  equation: 
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Equation  13 


where  t  is  ths  Incident  field  intensity  In  free  space 

i0  Is  the  reference  Inverse  distance  field  intensity 
T  is  the  transmission  distance 
D  is  the  distance  along  ray  path 

/(A)  Is  tha  functional  relationship  between  the  transmitter  dis¬ 
tance,  the  ray  distance,  and  the  radiation  angle 
A  is  tha  angle  between  ray  and  ground  plane  foran  inverse  distance 
field  intensity  of  300  millivolts  per  meter  at  one  kilometer 


Flgu  ■%  56  is  a  nomographic  solution  for  this  equation. 


(3>  Ground  reflection  losses  were  calculated  for  randomly  polarized 
sky-waves  for  various  angles  with  the  ground  for  three  types  of  earth, 
good  earth  U  •  10,  o  ■  I0'1  mhos/m),  poor  earth  («*4,  cr  «  lo“s  mhos/m), 
and  seawater  (t  »go,  tr  ■  5  mhos/m),  as  shown  in  Figures  57,  58,  and  59. 
Equal  amounts  of  energy  are  assumed  in  the  horizontal  ly  and  vertically 
polarized  fields.  These  reflection  losses  are  given  by  the  following 
equat I  on; 


L 


10  lot 


10 


Equation  14 


where  L  Is  the  earth's  reflection  loss  In  decibels 

RT  is  the  magnitude  of  the  vertical  polarization  reflection 
coefficient 

Rg  is  the  magnitude  cf  the  horizontal  polarization  reflection 
coefficient 

c.  Comparison  Between  Mea«uned  and  Calculated  values.  The  theo¬ 
retical  receiver  input  .voltages  based  on  free  space  anv.  ground  reflection 
losses  alone  were  calculated  for  the  various  probable  ray  angles  for  all 
the  circuits  with  sufficient  circuit  Information.  These  theoretical 
values  and  observed  yearly  peak  nighttime  receiver  input  voltages  agreed 
closely  for  -cst  of  the  circuits.  The  yearly  peak  receiver  input  voltage  io 
Hit!  ■•ig-itime  \alue  exceeded  ID  times  during  the  month  that  the  highest  val¬ 
ues  were  ooserved.  Taoic  5  shows  the  comparison  oetweencalculated  val jes  and 
values  observed  for  various  percentages  of  the  time.  The  statistical 
analysis  Indicates  that  the  calculated  voltage  based  on  f-ee  space  and 
ground  reflection  losses  alone  is  3.4  decibels  above  the  yearly  median 
value  of  the  composite  of  all. circuits  during  the  night  hours.  This  3.4 
decibels  value  can  bo  assumed  to  be  a  factor  to  obtain  median  field 
intensities  from  calculated  values  during  the  night  hours.  During  the 
day  hours,  the  additional  daytime  absorption  loss  must  be  subtracted. 
The  3.4  decibels  value  may  contain  additional  losses  due  to  magneto- 
ionic  splitting  and  layer  reflections.  For  many  circuits,  this  value 
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may  be  too  large.  These  circuits  have  a  single  predominate  sky-wave 
mode  throughout  the  entire  year.  Table  6  shows  comparisons  between 
calculated  and  measured  receiver  input  voltages. 

d.  Convergence  Gains.  In  comparing  theoretical  and  observed  values, 
theoretical  Ionosphere  convergence  gains  were  considered.  Figures  60 
through  63  show  convergence  gains  computed  for  reflection  of  electro¬ 
magnetic  waves  from  a  smooth  ionosphere  as  indicated  by  K.  A.  Norton.7 
The  observed  differences  between  the  calculated  and  measured  values 
indicated  that  the  convergence  gains  rould  generally  give  too  high 
received  voltages  and,  t*-us,  the  effects  of  convergence  gains  were  not 
included.  If  tho  magnitudes  along  any  ray  path  are  determined  by  areas 
only  a  few  Fresftel  zones  in  diameter,  convergence  gains  which  are  based 
on  reflections  over  a  wide  smooth  reflecting  area  are  probably  not 
applicable.  Scattering  and  wave  divergence  at  the  ionosphere  were  not 
computed. 

5.  Proposed  Prediction  Method  toi  Monthly  ..edian  Sky-wave  Field 
Intensities. 

a.  General.  The  proposed  prediction  method  is  designed  for  distances 
up  to  abou*  16.000  kilometers  for  short  wave  frequencies.  The  field 
intensity  values  obtained  are  for  tne  iu.ident  waves  from  the  ionosphere 
and  do  not  include  ground  reflected  components  at  the  receiving  location 
and  are  median  unabsorbed  field  intensities  less  ionospheric  absorption. 
An  unabsorbed  field  Intensity  is  that  expected  during  hours  of  darkness. 
These  values  given  in  this  report  are  for  a  reference  antenna  which 
radiates  a  signal  of  300  millivolts  per  meter  at  I  kilometer.  Unabsorbed 
field  intensities  for  any  antenna  can  be  obtained  from  those  based  on 
the  reference  antenna  by  accounting  for  the  difference  in  radiation 
between  the  actual  transmitting  antenna  and  the  reference  antenna.  If 
the  antenna  gains  are  calculated  with  reference  to  the  gain  of  a  short 
dipole  with  one  kilowatt  of  input  power  ( 300  millivolts  per  meter  at 
one  kilometer),  and  the  ai  tenna  radiated  powers  are  given  referred  to 
one  kilowatt,  the  unabsorbed  field  intensity  for  any  antenna  in  decibels 
is  the  transmitting  antenna  gain  plus  the  transmitter  power  added  to  the 
median  unabsorbed  field  inten&ity  of  the  reference  antenna.  Unabsorbed 
field  intensities  for  the  reference  antenna  arc  obtained  from  the  free 
space  distance  attenuation  along  the  ray  uaths,  ground  reflection  losses, 
and  an  experimental  factor  which  converts  these  theoretical  values  to 
median  values.  These  values  are  calculr;ed  for  angles  above  the  ground 
plane  which  are  determined  by  the  minimum  virtual  layer  heights,  the 
transmission  distance  along  the  earth's  surface  between  the  transmitter 
and  receiver,  and  the  number  of  reflections  between  the  ionosphere  and 
the  ground.  Sky-wave  field  i ntens  i  t  i  es  must  be  computed  in  slightly 
different  manner  for  short  and  long  distances.  The  method  for  distances 
short  enough  that  the  ionospheric  char --terist ics  at  the  path  mid-point 
determine  propagation  conditions  is  shn>  first  end  then  tne  method  for 
longer  distances  is  described. 


.  (  l ) 
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incident  down-coning  sky-wave  field  intensities  can  be  calculated  if 
the  transmitting  antenna  output  power  and  gain,  the  ray  angle  above  the 
ground  plane,  and  tne  dfstance  between  the  receiver  and  transmitter  arc 
known.  The  ray  anole  determines  the  antenna  radiation,  the  inverse  dis¬ 
tance  losses,  and  the  ground  reflection  losses  (for  routes  with  more 
than  one  ionosphere  reflection  point). 

1 2 )  A  graphical  procedure  Is  followed  requiring  various  charts  drawn 
on  the  same  rectangular  grids  astheCentral  Radio  Propagation  Laboratory 
monthly  predictions.  The  layer  height  is  found  first  and  then  using  the 
layer  height,  the  t  ransmlss Ion  d  istance,  and  the  number  of  hops;  the 
radiation  angle  is  determined  from  a  nomogram.  Average  F2  layer  Minimum 
virtual  heights  can  be  found  for  the  required  time  and  geographic  loca¬ 
tion  for  any  month  from  figures  23  through  34.  Figure  83  is  a  world  map 
and  Figure  84  is  a  great  circle  chart  used  in  determining  the  position 
of  the  great  circle  path  between  any  two  points  on  the  tarth. 

(3)  An  overlay  is  required  in  addition  to  these  charts  in  order  to 
determine  F2  layer  heights  for  any  hour  referred  to  any  reference  longi¬ 
tude.  A  transparency  is  placed  on  the  world  map  and  marks  are  placed 
over  the  position  of  the  transmitting  and  receiving  stations.  A  hori¬ 
zontal  line  isdrawn  over  the  equator  lineonthe  world  map  and  a  vertical 
is  drawn  at  the  longitude  for  which  the  hours  of  time  are  desired  to  be 
based.  If  it  is  desired  to  have  the  hours  in  Eastern  Standerd  Time,  the 
vertical  line  is  drawn  at  a  longitude  of  73  degrees  west.  This  overlay 
is  now  olaced  on  the  great  circle  chart  with  the  equatorial  lines  coin¬ 
ciding.  The  overlay  is  moved  horizontally  until  the  transmitting  and 
receivinq  locations  lie  on  the  same  great  circle  contour  or  at  propor¬ 
tional  distances  between  adjacent  qreat  circles.  Tne  great  circle  is 
sketched  in  and  the  path  midpoint  is  marked,  using  the  dash-dotted  scale 
of  the  great  circle  chart  as  the  distance  scale.  The  overlay  is  now 
ready  for  use  in  determining  the  layer  heights  for  any  month  for  any 
hour  referred  to  the  time  reference  for  the  particular  circuit.  This 
overlay  Is  placed  on  the  F2  layer  minimum  virtual  height  chart  for  the 
required  month.  The  equatorial  line  on  the  overlay  is  aligned  with  the 
equatorial  line  of  the  height  chart  and  the  vertical  reference  line  is 
placed  over  the  desired  hour  of  the  Say.  The  F2  minimum  virtual  layer 
height  is  read  at  the  path  midpoint. 

(4)  The  radiation  angle  for  any  assumed  number  of  hops  can  now  be 
determined  since  the  layer  heights  and  transmission  distances  are  known. 
Nomogram,  Figure  85,  is  used  in  this  determination.  Figure  e6  in  a 
nomogram  that  can  be  used  for  single  hop  modes.  These  nomograms  are 
based  upon  sharp  reflections  and  a  curved  earth.  The  slider  of  Figure 
85  is  olaced  ins  horizontal  position  at  the  layer  height  with  the  trans¬ 
mission  distance  scale  of  the  slider  set  at  the  left  vertical  reference 
line  of  the  nomogram  at  the  transmission  distance.  Thus,  the  radiation 
angle  can  be  read  directly  corresponding  to  any  number  of  hops  for  any 
height.  In  this  method,  all  E  layer  reflections  are  assumed  to  occur 
at  I  19  kilometers  while  the  F2  layer  reflections  occui  at  heights  based 
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cn  the  average  ■>*  f ccr  years*  »  mm. 

(5)  Knowing  the  radiatic.  angle  foi  »hc-  probable  mode,  it  is  now 
possible  to  find  tr,e  un*fcsorbed  field  intensity,  homonram,  rin.jro  -f, 
shows  the  free  spact.  ficid  *o.  n>  iitemi-  ■;  antenna.  This  nomogram  is 
entered  with  the  raciaticu  anjie  ana  ;!.,.  *  i-u  omission  distance  to  obtain 
the  free  space  field  for  any  my  pa.h.  7h<-  .mabsorhed  field  intensity 
for  the  reference  antenna  includes  odd  i  Mona!  losses  due  to  intermediate 
ground  reflections  and  a"  «»}  e«-!ircnft!  iac*ot  converting  the  calculate 
values  to  median  valu<=.  i.omog >.=  :;g,.reo  b"  through  53,  are  used  tc 
determine  the  ground  refection  losses  over  three  types  of  terrain,  good 
earth  (c  *  10,  <z  *  10-7  nhos/m  ,  r-‘o.  earth  !<:  =4,  <r  «  !0"3  mhos/m),  and 
seawater  (c  *80,  <7  ®  5  mhos/m),  for  ai  y  number  of  orouna  reflections  for 
any  radiation  angle.  ihe  strive.  :  or.  between  calculated  and 

median  values  was  shown  to  he  decitu-.  in  cur.rr.ary,  the  un- 

ebsorbed  field  intensities  for  toe  reference  r.« tem-a  are  obtained  by 
finding  the  mini  mum  vir‘ua.  layer  neign*s  from  a  nomogram,  finding  the 
radiation  angle  correspond  i  ns  to  any  number  of  hops  from  3  nomogram,  »n<i 
then  finding  the  free  spa«.e  field  intensity  from  another  nomogram.  The 
unabsorbed  field  intensi'y,  r-./.r*  r?  field  intensity,  for  the  refer¬ 
ence  antenna  is  equal  tc.  the  fret  sp.-^e  value  for  the  particular  trans¬ 
mission  distance  and  angle  eoove  the  ground,  less  the  ground  losses  at 
the  intermediate  ground  reflections  for  that  ladiation  angle  and  lose 
the  factor  of  ..!  decibels,  iteeerai'y,  one  iay  path  alone  determines 
the  sky-wave  field  intent  it\  for  most  distances  since  the  ionospheric 
conditions  limit  the  number  ■  '  r--;s:h|..  ray  paths  and  the  antenna  radia¬ 
tion  ii  generally  stronger  a  ong  one  ray  r.-.rn  in  addition  to  the  ground 
reflection  losses  ef'ectinn  toe  •*  '•■lgtii  ,f  multiple  hop  rays. 

c.  Daytime  FielP  Ihtins.tie'  ill  Ihe  daytime  incident  sky-wave 
field  intensities  ore  obtained  from  ui.<->'sorned  field  iricnsities  by 
subtracting  the  dayt  •  ITiC  flb  S><  option.  >■  o.-a.hica!  mt'hoo,  based  on  the 
semi-empirical  forr.ula  ;  the  statis*-cal  aoa'ysis  -s  used  to  calcuM** 
the  absorption.  it  is  required  to  know  the  solar  zenith  angle  at  the 
path  midpoint,  the  running  average  sunsrot  number,  the  angle  of  incidence 
at  a  100  kilometer  height,  t*'e  i>  wuer  of  ref  ections  between  the  g-tgund 
and  the  Ionosphere,  and  the  gy ro-f reuuuncv.  These  are  obtained  graphi¬ 
cally  from  the  standard  rectangular  qrids  mentioned  previously.  The 
angle  of  incidence  is  expressed  ir.  a  relationship  with  th&  rad'ation 
angle  determined  by  the  transmission  distance,  »he  layer  height,  and  me 
number  of  ionospheric  hops.  The  M.-aiing  average  sunspot  number  is 
predicted  while  the  wave  freauenc*  is  -tinii;  known. 

(2)  Tne  solar  zenith  angle  and  qy > u- f reouency  are  obtained  by  using 
an  overlay;  the  solar  zenith  a-'qle  cha-t.s,  f  .  lures  2  thro-igh  I?;  and  the 
gyro-frequency  chart,  figure  ml.  c  ■  overlay,  which  was  discussed  in 
proceeding  paragraphs,  is  placed  e  1  l  4j'»»  .'ooith  angle  chart  for  the 
required  month  with  the  equitt-’- :a ■  !  _•  molding  and  the  vertical 

reference  line  set  at  t*-e  dec.  i  rod  ic,i:  ■  he  solar  zenith  angle  is  read 
at  the  path  midpoint. 
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(3)  In  order  to  simplify  the  graphical  method,  a  new  function  Is 
defined.  This  function,  called  the  absorption  Index,  I,  Is  given  by 
the  following  formula, 

I  =  (COB  .  881  <P  )*'**  ( 1  ♦  •  0037  S  )  Equation  15 


where  S  Is  the  running  average  sunspot  number 
<p  Is  the  solar  zenith  angle 
I  Is  the  absorption  Index 

With  the  value  of  solar  zenith  angle  obtained  In  the  preceding  para¬ 
graph  and  the  predicted  running  average  sunspot  number,  the  Index  I 
can  be  obtained  from  the  nomogram  Figure  87,  The  Index  Is  read  at  the 
intersection  of  the  absorption  Index  scale  with  a  straightedge  placed 
between  the  solar  zenith  angle  and.  the  running  average  sunspot  number. 

14)  The  sky-wave  absorption  is  a  function  of  the  absorption  index, 
sum  of  the  wave  frequency  and  the  gyro-frequency, and  the  radiation  angle. 
The  distance  factor,  sec  <p,  Is  uniquely  determined  by  the  radiation  angle 
and  can  be  Implicitly  given  In  a  relationship  with  this  angle  above 
ground.  The  absorption  can  be  found  for  any  single  nop  transmission 
mode  from  the  absorption  Index  value,  the  radiation  angle,  sum  of  wave 
frequency,  and  the  gy ro-f requency  from  nomogram.  Figure  88.  This  pro- 
portjonal  nomogram  requires  two  settings  of  a  straightedge,  a  straight¬ 
edge  is  placed  between  the  absorption  index  valve  and  the  radiation 
angle  value  anda  mark  is  placed  onthecentral  index  line.  The  straight¬ 
edge  Is  then  placed  between  this  point  and  sum  of  the  wave  frequency  and 
the  gyro-f requency  and  the  absorption  is  read  on  the  absorption  scale. 
The  absorption  for  multi-hop  routes  Is  equal  to  the  absorption  for  one 
hop  routes  multiplied  by  the  number  of  hops. 


d.  Mode  Determination.  Die  extent  of  ionization  can  be  used  to 
eliminate  'mprobahle  modes.  It  has  beer,  shown  that  the  maximum  frequency 
usable  over  routes  of  various  lengths  is  proportional  to  the  secant  of 
the  angle  of  incidence  at  the  ionized  layers.  This  relationship  has 
been  incorporated  into  Nomogram,  Figure  89,  which  can  be  used  to  snow 
whether  or  not  a  layer  will  support  a  radio  wave  at  various  angles 
referred  to  the  ground  plane.  The  monthly  predictions  issued  by  the 
Central  Radio  Propagation  Laboratory  for  the  F2  4000  MUF  and  the  E  2000 
MUF  are  used  with  th*s  nomogram.  The  slider  of  the  nomogram  is  set  in 
a  horizontal  position  at  the  layer  height  with  eltner  the  F2  4000  MUF 
value  or  the  E  2000  MUF  value  (read  from  the  predictions  at  the  path 
midpoint)  on  the  slider  being  set  at  the  left  vertical  index  line.  The 
maximum  possible  radiation  angle  is  read  from  the  angle  curves  at  the 
wave  frequency  on  the  frequency  scale  of  the  slider.  The  value  of  mini¬ 
mum  F2  layer  virtual  height  used  is  that  found  in  the  preceding  para¬ 
graphs,  while  the  E  layer  neight  is  assumed  to  be  110  kilometers.  Thus, 
the  greatest  angles  for  E  layer  and  F2  layer  transmission  can  be  slmpiy 
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determined.  in  ordei  for  a  F2  layer  transmission  mode  to  be  possible 
during  the  daylight  hours,  *  he  waves  at  that  radiation  angle  must  be 
abie  to  penetrate  the  E  layer.  tms  nonograT,  can  also  be  used  to  find 
these  E  layer  penetration  frequencies.  1ft  slider  is  set  at  a  NO  kilo¬ 
meter  heigr.t  and  if  is  slid  horizontally  until  the  E  2 000  MUf  is  at  the 
left  vertical  index  line.  !he  t  layer  penetration  frequency  for  a  ^2 
layer  transmission  mode  is  read  on  slider  at  the  angle  for  the  F2  layer 
transmission  mode.  A  F2  layer  mode  must  penetrate  the  F.  layer  and  not 
exceed  the  F 2  'aver  maximum  usable  frequency.  Durinq  the  night  hours, 
it  is  only  necessary  to  find  the  maxiii-.-'"  radiation  angle  for  F2  layer 
modes . 

e.  Long  f'lstANCES.  rp.  path  lengths  exceeding  400C  kilometers,  there 
are  certain  modifications  necessary  in  tne  field  intensity  method.  The 
absorption  index  cannot  re-d  at  the  patn  midpoint.  An  average  ab¬ 
sorption  index  must  be  stu.  1 ne  solar  zenith  anqle  is  read  at  uniform 
inteivais  alony  tne  meat  .  ir>.le  path  between  points  2000  kilometers 
from  each  end  of  the  path.  The  absolution  index  is  found  sort  opjnu  i  ng 
to  these  solar  zenith  angles  and  averaged.  Tor  long  distances,  the 
Obliquity  factors  appear  "v  be  nearly  independent  of  the  number  of  nous. 
Thus,  the  absorption  can  be  solved  foi  any  angle  which  is  probable, 
inese  anqles  are  aqain  determined  I • om  somoqram.  Figure  8p,  which  shows 
radiation  angle  correspond i nu  •  o ary  transmission  distance,  layer  height, 
and  number  of  hops.  The  layer  heiqht  chosen  is  taken  to  he  the  average 
Of  the  heights  along  the  oreat  circle  pat*-  between  2000  kilometer  r  ‘ints 
from  the  ends  of  the  paths.  Oaths  wits  a.  lower  rjnber  of  hors  i  icwer 
radiation  angle)  are  favored,  due  *o  grtund  r.-flection  losses,  p''*!cu- 
lariyu/vet  poor  ground,  i.o.t  .tram,  lirjuie  89,  car  again  be  used  to  deter¬ 
mine  the  maximum  radia*ioi  anqle  and  to  determine  the  t  is.er  penetra¬ 
tion  angle  for  the  wave  'reqoen'y.  there  are  occasions  when  F;  layer 
propaqatior  is  improbable  a  ary  anqle  for  r  given  wave  (•  .er,  .  r.o 
critical  frequencies  are  re.fl  from  t*-  e",sel  ■■.•fro  *’•■*  --  '..'iP'ra- 
tory  predictions  for  the  reaui  red  "O-.tn  at  2000  x:,a  etei  ."i>  ts  from 
the  ends  of  the  path  to'  t  .e  -2  ‘Oyer  .  nd  122  ‘  iilonefer  points  from 
the  ends  of  the  path  for  the  l  layer.  h-  lower  value  at  the  two  path, 
ends  determines  the  maximum  usable  freot.*t  ies. 

f.  Examples,  several  ri'en*'  of  this  repnr*  were  used  to  snow  -cm- 
pa  r  i  sons  between  observed  .monthly  •  tdian  receiver  input  voltaqex  and 
calculated  values  uu<ii..c-u  by  tne  proposed  >  r  ••*»  r  *  *  n  •  inetnoo.  ih- 
results  are  aiven  on  rigurcs  ur.  •  ..ouu'-  u.s.  •.  ^mpie  ut  the  recuir.-q 
calculations  fur  the  proposed  prediction  netl.od  is  bi-.or.  in  Table  7. 
A  flow  chart  of  the  required  calculations  is  x~e.\~  i-  -ig-.re  94. 


6.  Conclusions 


a.  Absorption.  It  is  ne'ieved  »t. 
able  than  in  precedina  analyser  of  t*.r 
method  differs  f ron  cthe'r  "n  ... 

variations  in  ionospheric  layer  he  i  ; 
variation  between  seas.  is  account*  i 


'•r*.  “xi  menta  1  da»;.  was  avail- 
*.  vf«.  In  this  report,  the 
'*.*emt*  is  made  to  account  for 
h‘p.  'ne  anamol ooc  absorption 
f  x-  bv  layer  heigr-t  ja  r  i  at ;  ons. 
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The  secant  distance  factors  ara  computed  for  ionosphare  angles  of  inci¬ 
dence  at  a  100  kilometer  height,  roughly  corresponding  with  the  theoreti¬ 
cal  absorption  height.  The  gy ro-frequencies  are  also  computed  for  this 
100  kilometer  heJght.  The  analysis  of  the  measured  data  shows  that  the 
following  seml-empIrlcal  equation,  similar  In  form  to  Appleton's,  can  be 
used  to  account  for  daytime  absorption  at  frequencies  greater  than  about 
3  megacycles  per  second. 


A 


615.  S  n  see  <t 


[1  *  .0037  a] 

If  ♦  fa )S'9* 


[cos  .8auri0 


The  effect  of  solar  activity  on  absorption  is  based  on  the  calculations 
for  eight  circuits,  the  correlation  coefficients  ranging  from  .603  to 
•949.  The  correlation  coefficient  of  the  relationship  between  wave  fre¬ 
quency  and  absorption  is  .9424,  based  on  the  83  circuit  years  of  data. 

b.  Dav-to-Plv  variations.  Forty-three  circuit  years  of  mass  plots 
were  available  for  the  study  of  the  long  statistical  variations  of 
hourly  median  field  intensities.  The  difference  in  the  day-to-day  dis¬ 
tributions  for  a  period  of  a  month  is  quite  large  between  different 
months  and  circuits.  When  the  day-to-day  analysis  is  made  for  a  year, 
the  distributions  become  mere  nearly  alike,  the  variations  that  occur 
from  month-to-month  are  usually  explainable  by  seme  physical  cause,  such 
as  variation  of  ionospheric  support.  Vhen  the  data  is  analysed  over  a 
long  time  period,  both  the  day  and  night  voltages  have  distributions  that 
agree  at  least  roughly  to  the  Rayleigh  distribution. 

c.  absolute  Magnitudes.  Theoretical  nighttime  receiver  input  volt¬ 
ages  based  on  distance  and  around  reflection  losses  alone  agree  closely 
with  yearly  peak  nighttime  observed  voltages  for  most  of  the  circuits. 
Monthly  median  houriy  median  nighttime  voltages  were  obtained  from  the 
theoretical  peak  voltages  by  accounting  or  the  difference  between  peak 
and  meoian  hourly  median  values  based  uDOn  the  statistical  analysis  of 
the  day-to-day  variations.  Although  the  antenna  gains  could  not  be  cal¬ 
culated  with  greet  accuracy  for  some  of  tne  circuits  due  to  non-standard 
antenna  installations,  it  is  oelieved  that  some  of  the  discrepancies  are 
averaged  out  since  many  circuits  are  used  in  this  analysis.  No  arbitrary 
theoretical  factors  are  reqoirfed  in  computing  monthly  median  nighttime 
vo  I  tages. 

d.  Field  Intensitt  methc.  As  indicated  by  the  examples,  this  field 
intensity  calculation  method  is  reasonably  accurate  for  both  relative 
and  absolute  values.  Absorption  variations,  which  cannot  be  accounted 
for  by  variations  of  ihe  so'ar  zenith  angle  alone,  are  accounted  for  by 
variations  in  ionospheric  layer  heights  instead  of  presently  used 
arbitrary  seasonal  Tactors.  There  are  no  arbitrary  theoretical  factors 
in  this  method.  It  is  easy  and  fast  to  perform  since  all  operations  are 
graphical. 
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